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Pion String evolving in a thermal bath 
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By using the symmetry improved CJT effective formalism, we study a pion string of the 0(4) 
linear sigma model at finite temperature in chiral limit. In terms of the Kibble-Zurek mechanism we 
reconsider the production and evolution of the pion string in a thermal bath. Finally, we estimate 
the pion string density and its possible signal during the chiral phase transition. 
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I. INTRODUCTION 

Topological defects produced at phase transition are known to play an important role in early universe 
cosm ologyQBH- On the one hand, the early evolution of the universe undergoes a sequence of phase transi¬ 
tions and topological defects in these phase transitions may have observable consequences to the properties of the 
universe today. Such as the cosmic strings, they have been suggested as one possible source for the primordial density 
perturbations that give rise to the large-scale structure of the universe and the temperature fluctuations of the cosmic 
microwave background (CMB) radiationjlj], they could also contribute to structure formation or generate primordial 
magnetic fields which are coherent on cosmological scales 0- On the other hand, some predictive topological defects 
in the particle physics models, such as domain walls, which have problematic and cause a conflict with cosmology, 
can be ruled out[5|. 

The linear sigma model Q for the phenomenology of QCD has been proposed to describe the vacuum structure 
with incorporating chiral symmetry and its spontaneous breaking. The model can be used to describe a chiral phase 
transition in quantum chromodynamics (QCD) at finite temperature 0-El within the Cornwall-Jackiw-Tomboulis 
(CJT) formalismpjl. Similarly to the standard model of particle physics, traditionally, the linear sigma model does 
not give rise to topological defects which are stable in vacuum. However, if certain of the fields were constrained to 
vanish, it is still possible to construct a string-like configuration in the linear sigma model, the pion string, which 
would be topological defects [lil]. This kind of string is usually named as “embedded defects” [20] , which is not 
topological stable because any field configuration can be continuously deformed to the trivial vacuum. In order to 
cure this problem, a plasma stabilization mechanism had been proposed in Ref. [ED] [22f . In their works, they argued 
that the interaction of the pion fields with the charged plasma generates a correction to the effective potential and 
this correction reduces the vacuum manifold S 3 of the zero temperature theory to a lower dimensional submanifold 
S 1 , which makes the pion string stable. Moreover, this suggested stabilization mechanism for the pion string has been 
put on a firmer confirmation in recent works by Karouby and Brandenberger [23j [24j . So that wc have confidence to 
believe that pion strings are expected to be formed during the QCD phase transition in the early universe as well as 
in experiment of the heavy ion collisions. 

In order to study whether the pion string exists during the chiral phase transition, e.g. in the heavy ion collisions 
or in early universe. Firstly, we need to know how to estimate the production rate of pion strings in chiral phase 
transition, thereafter, we should know how to study the evolution and decay of these pion strings in a thermal bath 
when cooling from high temperature. Finally, we can reveal their possible trails in the experiment observations. 
So far the theoretical scenario that can be applied to study the formation of topological defects in rapid phase 
transition with global symmetry is the Kibble-Zurek mechanism [2§ - [28j . However, the CJT effective action in chiral 
limit presented in Refs. 1 71 - 1 171 | violates the Goldstone theorem and gives massive pions in the spontaneous symmetry 
breaking phase, besides, the temperature-dependent order parameter predicts a first order phase transition. All these 
not only disagree with the rigorous universality arguments, where the chiral phase transition can be of second order if 
the U(1)a symmetry is explicitly broken by instanton for Nf = 2 flavors of massless quarks [29j], but also conflict with 
recent results reported by Chiku and Hatsuda[30| [31) where by employing the optimal perturbation theory (OPT) 
and the real time formalism, they proved that the Goldstone theorem is always satisfied in any given order of the 
loop expansion in OPT in the 0(4) linear sigma model. Therefore, the direct use of the Kibble-Zurek mechanism to 
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the model has broken down. 

In order to restore the Goldstone theorem in chiral limit, Pilaftsis and Teresi recently developed a novel symmetry 
improved CJT formalism by consistently encoding global symmetries in loop-wise expansions or truncations of the 
CJT effective actionf32:]. Furthermore, this formalism has proven to avoid the existing problems and keep the features 
of the second order phase transition in the 0(4) linear sigma model[Hj. Thus, in this article we will apply the Kibble- 
Zurek mechanism to study the formation and evolution of the pion string embedded in a quark-gluon plasma at the 
LHC Pb-Pb collisions with a beam energy at ^/snn = 2760GeV by adopting the symmetry-improved CJT effective 
potential. 

The organization of this paper is as follows. In the next section we give a brief review of the pion string in the linear 
sigma model. In Sect. Ill, we calculate the effective potential by using the symmetry improved CJT formalism and 
obtain new equation of motion for the pion string. In Sect. IV, by using the Kibble-Zurek mechanism, the production 
and evolution of the pion string at finite temperature have been discussed. In Sect.V we discuss decays of pion strings 
for both the stable and unstable cases and we also address their possible observational consequence in experiment. 
Finally, Section VI summarizes our conclusions and discusses some open problems. 


II. THE PION STRING IN THE LINEAR SIGMA MODEL 

The Lagrangian density of the SU(2)r x SU(2)l symmetry linear sigma model has the form 


£ = - + d^if ■ d^if) - U(a,if), 


( 1 ) 


where the potential for the a and if is parameterized as 


U(a,if) = ^-(a 2 +n 2 ) + ^(a 2 +if 2 ) 2 


( 2 ) 


In our discussion of the pion string it proves convenient to define the new fields 


ip = 


a + in 


o 


V2 


( 3 ) 


and 


7T^~ = 


7T 1 ± in 2 

—JT' 


Then, the Lagrangian in Eq. 0 in terms of the p and fields is of the form 


( 4 ) 


£$ = + (d^TT+Xd^TT ) - XOXV + 7T+7T - ^ ) 


( 5 ) 


with x = A/6. The time-independent equations of motion read: 


V 2 ip = 2x{‘P*V + n + TT - y)(/? 


( 6 ) 


and 


„ f 2 

V 2 7t + = 2x{<p*<p + n + n — L ^-)n + . 


( 7 ) 


The global pion string extending linearly to the z-axis is obtained by solvin g th e equations of motion with an anti¬ 
symmetric vortex ansatz in the cylindrical coordinates in x-y plane, given by[. 


V = ^(r)e in0 , 

7T = * = = 0. 


(8a) 

(8b) 
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Here, the coordinates r and 9 are polar coordinates in x — y plane and the integer n is the winding number. In the 
following discussion, we will restrict ourselves to n = 1. Plugging the ansatz flSJ) into the above equations of motion, 
we get a second order differential equation 


with the boundary conditions 


P" + J ~ £ = X/ZiP* ~1)P, 


lim p(r) = 1, lim p(r) = 0. 

i—>oo r—>-0 


The energy per unit length of the string can be expressed as 


E = *fl 


drr 


p' 2 + ^ + $(p 2 - l) 2 


(9) 

( 10 ) 


( 11 ) 


where the first two terms come from the derivative of the field tp and the last term is from the potential of the a 
and 7f fields. Since the kinetic energy is logarithmically divergent, we need to introduce a cutoff R in Eg. (fill) which 
is the horizon size or the typical separation length between strings. Such a typical distance between strings should 
be determined by the initial string number density at the formation during the chiral phase transition based on the 
Kibble-Zurek mechanism. In general, R = O(fm). 

At tree level and zero temperature the parameters of the Lagrangian are fixed in a way that these masses agree with 
the observed value of pion masses and the most commonly accepted value for sigma mass. In the following numerical 
calculation, we take m a = 500 MeV and f n = 93 MeV as typical values. In the chiral limit, the coupling constant A 
is chosen to be A = 3 m 2 a / 


III. EFFECTIVE POTENTIAL COMPUTATION 


A convenient framework of studying chiral phase transition is the quantum thermal field theory. Within this 
framework, a resununation method introduced by Cornwall, Jackiw and Tomboulis (CJT) is an important and useful 
theoretical tooljTsj. However, in most of studies relevant to the truncated CJT formalism, there are two major 
drawbacks. One is that, in chiral limit, the CJT effective action violates the Goldstone theorem and gives massive 
pions in the spontaneous symmetry breaking phase. The other is that the temperature-dependent order parameter 
predicts a first order phase transition. This prediction, of course, disagrees with the rigorous universality arguments, 
where the chiral phase transition can be of second order if the U(1)a symmetry is explicitly broken by instanton for 
Nf = 2 flavors of massless quarks [29[. While the former conclusion comes into conflict with recent results reported by 
Chiku and Hatsuda[30| [3l|] , where by employing the optimal perturbation theory (OPT) and the real time formalism, 
they proved that the Goldstone theorem is always satisfied in any given order of the loop expansion in OPT in the 
0(4) linear sigma model. Therefore, in following discussions, we take the symmetry-improved CJT effective as a 
prototype framework. 

In the case of the linear sigma model, the CJT effective potential at finite temperature will appear asjl8| |32] [33f 


V{<l>,G a ,G„) = 


U{4>) + l -J\nG-\frk)+*-j\n G-\fr k) 

J [D- 1 (frk)G l ,(<frk)-i\ +\J [D~\(l>-,k)G*{hk)-l\ 


( 12 ) 


+V2(<t>, G a , G n ). 

Here, after shifting the sigma field as a —> a + (f>, the classical potential U((f> ) has the form 

u{<)>) = 

and the last term 14(<^j G a , G T ) denotes the contribution from the “ 00 ” diagrams, which is equivalent to the Hartree 
approximation, explicitly, 


(13) 


I'M'A Gu , Gjr) — 3 


A_ 

24 
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A_ 

'24 


GAfok) 




+1 4 


1 2 


G^k) 


(14) 
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Here and there, we use the imaginary-time formalism to compute quantities at nonzero temperature, our notation is 

/3 ?/ i] s l l> ’ 

where /3 is the inverse temperature, /? = jrj-f-, and as usual Boltzmann’s constant is taken as ks = 1, and = 2itnT, 
n = 0, ±1, ±2, ±3, • • •. For simplicity we have introduced a subscript /3 to denote integration and summation over 
the Matsubara frequency sums. Also, we denote D a and as the tree-level sigma and pion propagators, while G„ 
and D n are treated as the full propagators. By minimizing the effective potential with respect to full propagators, we 
obtain a set of nonlinear gap equations for the sigma and pion thermal effective masses 

Ml = m 2 + ±F{M a ) + ^(M w ), 

Ml = m 2 + ^4> 2 + ^F{M a ) + ^ F(M n ) 

DO 0 

Here we have used a shorthand notation and introduced the function 

F ( M ) = [ 2 r l „ f2 - (16) 

J p u 2 +k 2 + M 2 

In order to recover the Goldstone theorem in the truncated CJT effective action, in analogy to the standard One- 
Particle-Irreducible effective action, a symmetry improved CJT formalism has been proposed in Ref. [32| through the 
introduction of a constraint, which can be consistently implemented by redefining the truncated CJT effective action 
with a Lagrange multiplier field. It has the form 


(15a) 

(15b) 


ml = 0. (17) 

Substituting this constraint to the above gap equations urn in the symmetry breaking phase, the constrain (ED 
implies the massless pions are naturally realized whenever the order parameter is not vanished. On the contrary, in 
the symmetric phase, we have <j> = 0 and the constraint is automatically satisfied. As a consequence, the two mass-gap 
equations become degenerate, the particles have the same mass with M a = M„ = M, and 

M 2 = m 2 + \F(M). (18) 

In order to get more insight into the nature of the phase transition and verify that the transition is of the second 
order, We solve the system of mass-gap Eqs. ED with the constraint ED using a numerical method based on the 
Newton-Raphson method of solving nonlinear equations. In this way, we are able to determine the effective masses 
M a and M v and the order parameter cf> as functions of temperature T. 

Using the definition in the framework of the symmetry improved CJT formalism, we can rewrite the finite tem¬ 
perature effective potential 14 // as a function of t/j in the thermal Hartree approximation by extending the mass-gap 
equations (fl5l) from (f> i/; while 


Ml(if>) = - dFe//W , 
i\) dip ’ 

with 

\ K \ \ 

Mlty) =m 2 + —F{M a ) + — F(M„) + -ip 2 . 

ODD 

Furthermore, in chiral limit, in the symmetry breaking phase, M'p = 0, then 

Ml(i/>) = ——(j) 2 H— ip 2 . 


(19) 


( 20 ) 


( 21 ) 


The order parameter <p as functions of temperature T can analytically expressed in terms of the T c as pointed out in 
Ref, frill as 


V = l -(T 2 c -T 2 ) = fl{ 1-^J; 


( 22 ) 
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and the critical temperature T c = v^/V — 131.5MeV. Moreover, the effective mass of sigma meson in the symmetry 
breaking phase can also described as a function of the T c 


Mi 




(23) 


Now the Goldstone boson masslessness condition can be naturally implemented in the symmetry breaking phase, and 
the integration of the effective potential V e // with respect to the field ?/> really create a typical potential possessing 
the second order phase transition. 

Integrating the equation (1191) , we can determine the symmetry-improved effective potential 14// as a function of i[> 
in Hartree approximation as the form of 


Vefffy) = - </> 2 ) 2 , (24) 

or more specific in term of the tp field as 

Veffb) = x(<P 2 - y) 2 . (25) 

If we approximately treat a hot thermal medium as a uniform thermal bath with pion strings embedded in it, a new 
set of coupled equations of motion for the pion string could be derived by simply replacing the relevant classical 
potential with the above appropriately modified thermal effective potential 14//. Accordingly, in term of the ip field, 
the equation of motion of ip field in Eq.((6|) can be described as 

V 2 p = 2x(?V - yte- (26) 

From equations (1251) and (1261) . we can naturally recover the classical potential flSJ) and the equation of motion (O for 
the pion string at zero temperature. 

As mentioned in our previous study [331], the constraint proposed in Eq. lfTTl) could be obtained without any assump¬ 
tion in the large -N approximation in the framework of the usual CJT formalism, thus the main results presented above 
can also be applied to the case of the larg e-iV approximation directly only that the critical temperature T c = is 

replaced by T' c = — 1 6 1 MeV ITT} |ll[ and this critical temperature T' c is closer to the Lattice QCD prediction [34l— 

1371 ]. However, in this work, without any loss of generality, we prefer to treat the critical temperature T c as an input 
parameter in order to give out a more realistic analysis of the production and evolution of the pion string in high en¬ 
ergy heavy ion collisions, e.g. for LHC energies. The reason for this modification is to be as close to the experimental 
situation as possible. 

We first numerically investigate pion string solutions at finite temperature. In Fig[T| we plot the profile functions of 
the pion string for different temperatures. It is shown that all the fields are moving towards to the trivial values with 
increasing temperature. When T is larger than some critical temperature T c , there exist only the trivial solutions for 
the equation of motion in Eg. (1261) . and the pion string is melted away. Correspondingly, the trivial solution indicates 
the restoration of chiral symmetry in full space. 


IV. THE KIBBLE-ZUREK MECHANISM 

In order to study whether the pion string exists during the chiral phase transition, we have to estimate how many 
pion strings are created in the phase transition by the Kibble-Zurek mechanism |25 m 28| . As described by the Kibble- 
Zurek mechanism, the pion strings can be expected to form in a second order phase transition even it is a perfectly 
homogeneous fast transition, if the transition proceeds faster than the order parameter of the breaking symmetry 
phase is able to relax. In such a non-equilibrium transition, because of fluctuations of the order parameter, the new 
symmetry-broken-phase starts to form simultaneously and independently in many parts of the system. Subsequently, 
these independent regions grow together to form the new symmetry-broken phase randomly. At the boundaries, 
where different causally disconnected regions meet, the order parameter does not necessarily match and a domain 
structure is formed. Such a random domain structure reduces to a network of strings if the broken symmetry is the 
U(l) symmetry. 

Facilitated with the Kibble-Zurek mechanism, we consider a system cooling from an initial temperature T l through 
the critical temperature T c by the change of a control parameter e. For quantitative estimation, we take a quark-gluon 
plasma created at the LHC Pb-Pb collisions with a beam energy at y/S nn = 2760GeV as a thermal bath. As a result 
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FIG. 1: Numerical solution of the pion string in the case of chiral limit with n= 1. (a) The profile functions of pion strings as 
function of the radius r at various temperatures, where the solid curve is for zero temperature, the dashed curve is for T = 0.5 T c 
and the dash-dotted curve is for T = 0.977 . (b) The energy densities as a function of the temperature T. 


of ultrarelativistic collision between two heavy ions, a fireball is expected to form and then it will rapidly thermalize 
due to the fact that there exists manifestly partons with very large cross section for gluon scattering. The enormous 
amount of energy density deposited in the fireball results in large pressure gradients from the central to the peripheral 
region of the fireball that drives the expansion of the fireball. This expansion leads to cooling of the fireball and the 
hydrodynamic model gives that the temperature of the fireball T evolves with the time t as [38N42l ] 

T(f)ocfT (27) 

Furthermore, with the expansion of the fireball, the interparticle distance grows with time (or temperature), the 
particles cease to interact after sometime and free stream to the detector. The surface of last scattering is called as 
the freeze-out surface. Since scattering could be both elastic (where particle identities do not change) and inelastic 
(where particle identities change), it is possible to have two distinct freeze-outs, namely, chemical freeze-out, where 
inelastic collisions cease, and thermal/kinetic freeze-out where elastic collisions cease and the particle mean free 
path becomes higher than the system size, which forbids the elastic collision of the constituents in the system (43). 
In experiment, the chemical freeze-out surface is determined by analysing the measured hadron yields, while the 
kinetic freeze-out surface can be determined by studying the data of transverse momentum distribution of produced 
particles. In Ref.[44j], the chemical and kinetic freeze-out scenarios in central heavy ion collisions from the lower 
AGS energies to the largest LHC energies are investigated in detail. Based on their studies, the extracted chemical 
and kinetic freeze-out temperature corresponding to LHC energies for zero chemical potential are T c h — 150MeV and 
Tki n — 90MeV, respectively. The radius which is reflective of volume of the fireball can also extracted as r c h ~ 10.87fm 
at T c h ~ 150MeV. Correspondingly, the volume of the fireball at the moment is about V (T c f, ) ~ 5380fm 3 and the 
subsequent evolution of the volume is given by 

V(t) oc t. (28) 

Employed with the resulted thermal parameters, temperature T c h , we can now fix the critical temperature T c . It 
was realized that there are very high densities in the system when the temperature is close to the critical temperature 
T c , but as long as the temperature is below the critical temperature, the quark-gluon plasma created at heavy ion 
collisions experiences a hadronization, and the rapid fall-off of density lead to rapid equilibration and chemical freeze- 
out only a few MeV below the transition temperature (45| . The current understanding of the connection between 
the QCD phase diagram and chemical freeze-out can be revealed through the comparison of lattice calculations with 
results from the hadronic gas models. Therefore, from the recent work of [3f|, the precise critical temperature is 
T c ~ 157MeV for zero chemical potential. 

In the vicinity of T c , a continuous second-order phase transition can be characterized by the divergence of both the 
equilibrium correlation length £ 



( 29 ) 
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and equilibrium relaxation time r 

r{e) = ^7, (30) 

as a function of the distance to the critical point T c . Where £0 is the zero temperature-limiting value of the 
temperature-dependent coherence length £(T), and in this work we take £0 = 1 /m a — 0.4fm. But for To, since 
the experimental data obtained in the RHIC and LHC experiments favor a very short thermalization/equilibration 
time, we take ro ~ 0.06fm so as to gain a reasonable value for the critical time t c at T c ~ 157MeV for zero chemical 
potential. From Eg. (1251) . the correlation length £(T) is given by £(T) = (^A^ 2 ) -1 / 2 , where <j> is the global minimum 
of the effective potential in the homogeneous case [361 ]. Consequently, it is convenient to define our reduced distance 
parameter as 


T 2 - T 2 


£ 


I / 


= 1 - 


T 2 

T 2 ’ 


(31) 


and we can get the correlation length critical exponent v = 1/2. While the dynamic critical exponent z is equal to 
2 in a Ginzburg-Landau system at a second order phase transition. In terms of the reduced distance parameter e, 
the system initially prepared in the high-symmetry phase (e < 0) is forced to face a spontaneous symmetry breaking 
scenario as the critical point is crossed towards the degenerate vacuum manifold (e > 0). 

The reduced parameter is characterized by the quench time tq , 


T Q 



(1 

dT 


dt 


and varies linearly in time <; according to 


e = 


V 


(32) 


(33) 


in ? £ [— tq,tq], the critical point being reached at c = 0 and the system time t is given by t = t c + ?. Far away 
from the critical point T T c or T -C T c , the equilibrium relaxation time is very small with respect to the time 
remaining until reaching the critical point following the quench (1331) . and the dynamics is essentially adiabatic. In the 
opposite limit, in the close neighborhood of e = 0, the dynamics is approximately frozen due to the divergence of the 
equilibrium relaxation time (critical slowing down). The system is then unable to adjust to the externally imposed 
change of the reduced control parameter e, and the order parameter of the system ceases to follow the equilibrium 
expectation value and enters an impulse stage within the time interval [— t z ,t z \. This is often referred to as the 
adiabatic-impulse approximation, which is believed to capture the essence of the nonequilibrium dynamics involved in 
the crossing of the phase transition at a finite rate. The inability of the order parameter to keep up with the change 
imposed from the outside is the essence of the freeze-out. 

The above time boundary between the adiabatic and frozen stages t z can be estimated as follows. When <7 > 0 or 
the temperature is below T c , the order parameter coherence spreads out with the velocity 


c(T)~^ = ^ei (34) 

T T 0 

The freeze-out of pion strings are expected to occur at the Zurek time t z when the causally disconnected regions 
have grown together and the system is to establish coherence in the whole volume. In the same time, the global 
thermodynamic equilibrium of the system may be realized and the order parameter will follow its equilibrium value 
in whole space. Therefore, this moment is also corresponding to the end of the hadronization process in heavy ion 
collisions and the freeze-out temperature T z at the Zurek time t z in the Kibble-Zurek mechanism can be directly 
identified as the chemical freeze-out temperature T c h due to the fact that the very hadronization process leads to the 
equilibrium distributions of hadrons [4Tj. 

At the freeze-out temperature T z (T z = T c h), the causal horizon has increased to the distance 


Zh(s») = C ( T )* = (35) 

Jo r o 

which has to be equal to the coherence length £(^ z ) as shown in FigJ5] This condition leads to a relationship between 
the freeze-out temperature T, and the freeze-out time <; z as follows: 


2 — 



t q' 


e. 


(36) 
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FIG. 2: With decrease of temperature T < T c , the causal horizon will move out to a distance £h, and it will be subsequently 
equal to the coherence length £(? z ) when the equilibrium phase is achieved at the freeze-out point T z ~ 150MeV. 


Then, it is easy to obtain 


t z =t c + <; z =t c + \Jtqtq ~ 9.60fm, (37) 

when temperature is at T z ~ 150MeV. Here, t c = 3tq/ (^(T c /T z ) 3 — 1 ^ ~ 8.37fm, which is consistent with the 

calculation of the hydrodynamic models used to interpret the RHIC and LHC results^. 

The Kibble-Zurek mechanism sets the average size of the domains by the value of the equilibrium correlation length 

at £ e , 


6 = = £o( — )* - 1.81fm. 

to 


(38) 


Moreover, it is worth to point out that the equilibrium correlation length at e z should be taken as the horizon size or 
the typical separation length between strings R. In general, the correlation length £ 2 is often recast as an estimate 
for the resulting density of topological defects, 


1 


n = 




(39) 


where the parameter re is introduced to reduce the overestimation of real density of defects in experiment and k ~ 
1 — 100 depends on the details of the specific model. In this study, the experimental measurements in heavy-ion 
collisions will give a constraint on re. n is defined as string length per unit volume. We can estimate n ~ 0.3fm/fm 3 
for re = 1 and n ~ 0.003fm/fm 3 for re = 100. This is the main prediction of the Kibble-Zurek mechanism, and the 
result presented here can be generalized to other Ginzburg-Landau systems having a second-order phase transition. 


V. THE DECAY OF THE PION STRING 

Note that in our situation, closed string loops give dominant contributions to the total energy of the strings. This 
is because at the freezing out time t z , the string evolution is still ruled out by the frictional force by the surrounding 
matter so that the free motion of the string is not realized. In addition, the expansion of the system is too rapid that 
the initial Brownian string distribution will be conserved. Then the initial structure of the string network partially 
remain and the spatial trajectories of strings are very complicated. Thus we take the initial pion string network as 
that of the Brownian one and regard that the distribution of these loops does not change with time t <t z so long as 
they could survive in such a ruinous environment, the size of the loops are conformally stretched during the expansion 
of fireball and a simple scaling can be realized. 
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By using the scale invariance of Brownian string described by Vachaspati and Vilenkin [47}, the initial distribution 
of number density of the pion string with the length between l and l + dl at the chemical freeze-out temperature T c h 
can be described as 


dn(l) = h~Ul, 


(40) 


where the parameter k is not a free parameter and it is easy to establish the relationship between the above parameter 
k in equation (1391) with the parameter k here. As we already know n is the string length per unit volume and the 
volume of the fireball at T z is V (T,) ~ 5380fm 3 , then the total length of the pion string in the fireball at the freeze-out 
time t z is set by Lt — n x V(T Z ). On the other hand, as long as we obtain the distribution of number density of the 
pion string n(Z), the total length of the pion string in the fireball at the freeze-out time t z is given by 


L t = V(T Z ) / ldn(l ) = kV{T t )£ 


r*dl = 2kV{T z )Z z *{l a * -l b *). 


(41) 


Here, l a and l b are the minimum and maximum lengths of the pion strings, respectively. Note that the string width 0 
gives a minimum length of the pion string at time t z , l a — 2ir£ z — 11.37fm, while the longest string is also constrained 
by the volume of the system with lb — 27rr c / l ~ 68.26fm. After a simple calculation, we arrive at 

k = y ~r yyr - • (42) 

2/c£j(Z„ T -l b 3 ) 

Then if n = 1, k is about 2.12, whereas k is about 0.02 for k = 100. 

After integration of the equation flUl) . we have the total number of pion strings: 

N(T Z ) = V(T Z ) [ ' dn(l) = 2kV ^ z) (0 1 - Z"*) ~ 58 k. (43) 

Jia 30 V ' 

Then the total number of pion strings created at the LHC Pb-Pb collisions with a beam energy at ^/sjvat = 2760GeV 
could vary between N ~ 116 and N ~ 1. 

Unfortunately, in order to detect these produced pion strings in the heavy ion collisions experiments, we need 
to deliver a comprehensive analysis on the process of the pion string decay after the time t z . For simplicity and 
convenience, we separate our subsequent studies into two cases which are corresponding to two limits of the stabilities 
of pion strings. The first case is that the produced pion strings are suggested to be absolutely stable during the time 
from the chemical freeze-out time t z to the kinetic/thermal freeze-out time tf according to the plasma stabilization 
mechanism proposed in Refs. j2ll - l24l |. But, as soon as the time is across tf, the particles in hadron resonance gas 
depart from each other so fast that the collision processes become ineffective, thus pion strings become unstable and 
decay suddenly. On the contrary, there exists another possibility: the plasma stabilization mechanism does not work 
very well in the case of the heavy ion collisions experiment, and the produced pion strings during the chiral phase 
transition will decay immediately after their births when t > t z . We will point out later that the second case is closer 
to the experimental situation. 


A. The stable case 

In this case, after the freeze-out temperature T z , the pion string will stop forming and scale with the expansion 
of the system or the fire ball created in the heavy ion collisions until to the kinetic freeze-out time tf. Because the 
volume of the system obeys the following law V(t) ~ t , then we make an assumption that the correlation length of 
the pion string will rescale as 

&=&(5r)*> (44) 

until to the finial stage of the evolution of the fire ball. Based on the chemical and kinetic freeze-out scenarios in 
Ref. [44j . the kinetic freeze-out temperature will occurs at T/ ~ 90MeV, so that we have tf ~ 44.43fm. Now the 
correlation length of the pion string achieves at 0 ~ 3.02fm, which gives out the minimum length of the pion string 
at time fj as Lj = 2-7T0 ~ 18.97fm. 

By making the assumption that pion strings are absolutely stable in a plasma and can survive after the decoupling 
time t f , the evolution and decay of the formed pion string has been extensively investigated in Ref (27j . The qualitative 
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results obtained in that work are applied in this work too in the case of slight modification of the parameters over 
there. In order to avoid the simple repetition, in this section, we only bring forward some comments on the previous 
work presented in Ref 27 . 

The first comment: the thermal or kinetic freeze-out is the stage in the evolution of matter when the hadrons 
practically stop to interact, at that time, the transverse momentum distributions of the particles get fixed and do 
not change thereafter. If pion strings can survive after the decoupling time tf , neutral pions and sigma (charged 
pions) particles emitted from pion strings do not get chance to be thcrmalized. These unthcrmalized pions will give 
a reasonable contribution to the px spectra of 7r and induce a visible fluctuation/distortion in the low pr-momentum 
regimes. In fact, the experimental data show the hadron spectra in the soft region, pt < 2GeV, have thermal character, 
although the original thermal distributions are modified by the collective transverse flow and decays of resonances. 
So that if pion strings are really formed during the QCD chiral phase transition, they should decay completely before 
tf. Otherwise, they could be already observed in experiment. 

The second comment: based on the chemical and kinetic freeze-out scenarios in Ref. [44|, the chemical freeze-out 
time t z is the moment when the inelastic collisions stop, so that particle identities do not change and the temperature 
of the freeze-out is inferred from the studies of the ratios of hadron multiplicities. Moreover, when t > t z , the hadron 
resonance gas is taken as a gas of noninteracting hadrons and resonances which is a good low T c h approximation to 
QCD thermodynamics. Any additional pions burst from pion strings when T < T c h would pose a serious change 
to the standard chemical freeze-out models which have been quite successful in describing the particle multiplicities 
measured in heavy-ion collision experiments. Therefore, it is more reasonable to believe that pion strings will decay 
immediately soon after they are produced in heavy ion collisions. 

The third comment: the plasma stabilization mechanism proposed in Refs. [2lU24j does work only when the coupling 
constant A is extremely small. For a realistic value of A in this work, the temperature condition for the stability of 
the pion string is 


T>T d = ~ 2335MeV. (45) 

e 

Td is even much larger than the T c , so it is natural that pion strings produced in LHC experiment are unstable and 
they will decay quickly after giving birth. 


B. The unstable case 

In the Kibble-Zurek mechanism, the fast second-order chiral phase transition of the fireball leads to the formation 
of pion strings at freeze-out time t z . In analogy to a cosmic string created in cosmology [lj, the pion strings formed in 
nucleon-nucleon collisions may be imagined as sigma - neutral pion pairs confined in a string-like tube. In the next 
step, the strings decay/fragment forming directly hadrons (the sigma particles and neutral pions), and the produced 
hadrons are modeled as smaller pieces of the original string. 

As mentioned above, the pion string is not topologically stable, since any field configuration can be continuously 
deformed to the vacuum. It can be shown by numerical analysis that the pion string is only stable for very small 
values of A. This implies that the pion string is unstable in the parameter space allowed experimentally[l9j]. In the 
early universe and in heavy-ion collisions, pion strings are expected to be produced and subsequently decay. Their 
lifetime can be estimated by considering their interactions with the surrounding plasma. Based on a naive dimensional 
analysis, their lifetime r pa should be proportional to the inverse of the their masses/energies at specific temperature. 

From the right panel in Fig[l] the energy per unit length of the pion string can be parameterized as 

E(T) = Eo{l - ^J), (46) 

with Eq is the energy of the pion string at vacuum, and Eq ~ 231.38MeV/fm. At T = T Z , the energy per unit length 
of the pion string is about 20.17MeV/fm. Since a minimum length of the pion string at time t z is l a — 27t£ 2 ~ 11.37fm, 
then the mass of the pion string M(T ) at T = T Z at least has 

M(T Z ) = 2t t£ z E(T) ~ 229.3MeV. (47) 

Correspondingly, the lifetime r ps of the pion string at moment is less than 0.86fm. Thus we can see the pion string 
indeed decay immediately soon after they are produced in heavy ion collisions at t = t z . 

At chemical freeze-out, the inelastic interactions among the produced particle stop. However, these particles can 
still interact elastically which could affect their momentum distributions. Then all sigmas and pions burst from pion 
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strings will be thermalized quickly in a plasma. Furthermore, since the thermal character in the hadron spectra in 
experimental data can not be modified by the collective transverse flow and decays of resonances. It is extremely 
difficult to recognize and isolate sigmas and pions decayed by pion strings from the background thermal productions 
by using the efficient method of the transverse momentum distributions or the HBT interferometry (more recently, 
the name femtoscopy is used) technique (48}. 

Even though the produced pion string is unstable, there still have possible signals of the pion string produced because 
the string configuration violates the isospin symmetry. Since the generic expectation is that pions of different isospins 
would be produced in equal abundances because of isospin symmetry, however, the formation and decay of pion strings 
would produce abundant neutral pions, which would lead to a large deviation in the charge-neutral correlation from 
expectations based on the generic pion production mechanism. Similarly to the DCC (disoriented chiral condensate) 
phenomenon [49j - f5l] | . pion strings are also expected to emit pions coherently from the collision volume, which may 
results in large fluctuations (sometime called as DCC-like fluctuations) in the neutral pion fraction, /, defined as 
/ = no/(no + n c ), where no and n c are the multiplicities of the neutral and charged pions, respectively. 

Experimentally, the existence of DCC phenomenon was previously investigated in heavy-ion collisions at the CERN 
Super Proton Synchrotron at Jsnn = 17.3GeV fihil and at Tevatron in p + p collisions by the MiniMax col¬ 
laboration at y/sNN = 1.8TeV[56| 57].In both cases, the possibility of large-sized DCC domain formation has been 


excluded by the measurements. Recently, event-by event fluctuations of the multiplicities of inclusive charg ed p ar- 
ticles and photons at forward rapidity in Au + Au collisions at y /snn = 200GeV have been studied in Ref. [58j] , in 
their studies, it is shown that for all pions, the observational deviation is found to be less than 1% from the generic 
expectations based on the generic production mechanism of pions owing to isospin symmetry.This means that if pion 
strings or DCC domain are really produced in the heavy-ion collisions, the number of the produced pions should be 
less than 1% of the total pions in fireball. Under this constrain, at the LHC Pb-Pb collisions with a beam energy 
at yj snn = 2.76TeV[5g|, the numbers of sigmas and neutral pions produced from pion strings are expected to be 
N n o ~ N a ~ 25 — 7 if we ignore the effect of the DCC phenomenon or other effects. Then the total number of pion 
strings created at the LHC Pb-Pb collisions are about N ~ 4 and IV ~ 1, which give k £ [0.02, 0.07]. Therefore, the 
production rate of topological defects seems much lower than that of system in condensed matter physics jgO]- 


VI. SUMMARY AND DISCUSSION 

We have investigated the production and evolution of the pion strings in a thermal bath within the framework 
of the 0(4) linear sigma model in chiral limit by adopting the symmetry-improved CJT effective formalism. It has 
been shown that the formalism is keeping the Goldstone theorem and satisfying a second-order phase transition, this 
makes it possible to apply the Kibble-Zurek mechanism to study the formation and evolution of the pion string in a 
chiral phase transition in early universe or heavy ion collision. Our results indicate that pion strings are expected to 
be produced in a second-order chiral phase transition, and then quickly decay into neutral pions and sigmas after the 
chemical freeze out time t z — 9.6fm. Following the Kibble-Zurek mechanism, we can roughly estimate the number of 
pion strings are to be 1 — 4, and their subsequent decays will produce about 7 — 25 neutral pions and sigmas. Because 
the string configuration violates the isospin symmetry, the neutral pions are more likely to be produced during the 
chiral phase transition when comparing to the charged pions. This can be taken as a possible signal of the pion string. 

However, the possible signal observed in this work is very similar to the experimentally observable anomalies caused 
by the DCC decay. A further effort needs to be addressed in future in order to clarify the differences from the pion 
string decay and the DCC decay. Another problem deserved to be noted is that chiral phase transition is actually a 
crossover rather than a phase transition, then the estimation of the pion string number density is not so straightforward 
by directly using the Kibble-Zurek mechanism, and simulation of the production and evolution of pion strings in early 
universe or heavy ion collision remains as an important future problem. These two problems are certainly out of the 
scope of our current topic and we prefer leaving for future study. 

Before concluding, let us give some notes about the existence of topologically stable non-Abelian global strings in 
the U(N) linear sigma model j6li - l63j . If the axial anomaly is absent, there exist stable Abelian axial string (or rj 
string[19f ) winding around the spontaneously broken U(1)a and non-Abelian axial string winding around both the 
U(1)a and non-Abelian SU(N) chiral symmetries. Therefore, If the U(1)a symmetry is exact above the QCD chiral 
phase transition, then, as long as the effects of instantons can be neglected, the rj string and non-Abelian axial string 
can form during the chiral phase transition of QCD, and moreover, they are topologically stable. This is unlike the 
pion string, these produced Abelian and non-Abelian axial strings can survive until to the kinetic freeze-out time tf, 
when t > tf, the particles depart from each other very fast and the collision processes become ineffective. Thereafter, 
large pressure gradients from the central to the peripheral region of fireball suddenly switch off and these stable 
strings start to shrink to their minimum lengths then decay quickly. Moreover, as the density of matter and/or 
the temperature decreases, it is expected that the instanton effects will rapidly appear and the axial anomaly term 
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has to be included in the model. Therefore, metastable domain walls are present and Abelian axial strings must be 
attached by N domain walls, forming domain wall junctions[s] j64|. Thenceforth, Instable domain wall junctions will 
decay sooner or later through quantum and thermal effects. It is interesting to investigate whether domain walls and 
strings are stable or not until the kinetic freeze-out time tf since they seem to live much longer than pion strings. 
And we believe any kind of topological solitons created in the heavy ion collisions will be interested both for the 
theoretical and experimental studies. Because these formed stable solitons can be directly detected and confirmed by 
the experimental observation if they can live on until the kinetic freeze-out time tf. Eventually, work in this direction 
is under progress. 
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